Exclusive photoproduction at the Tevatron and LHC within the dipole picture 



OO 

o 
o 

(N 



(N 
> 



- '—I 

X 



L. Motykai'2 and G. Watt^ 

'II Institute for Theoretical Physics, University of Hamburg, 22761 Hamburg, Germany 
^Institute of Physics, Jagellonian University, 30-059 Krakow, Poland 
^Department of Physics & Astronomy, University College London, WCIE 6BT, UK 

We provide predictions for the rapidity distributions of exclusive photoproduced J/t/j and T 
mesons, and Z'^ bosons, at the Tevatron and LHC. We use the equivalent-photon approximation 
with the photon-proton cross sections given by the impact parameter dependent dipole saturation 
model, which has already been shown to give a good description of a wide variety of HERA data. 
We derive the quark-antiquark light-cone wave functions of timelike neutral electroweak bosons. An 
essential difference is pointed out between the amplitude for timelike heavy boson photoproduction 
and the amplitude for deeply virtual Compton scattering. 



I. INTRODUCTION 



Exclusive diffractive vector meson production, V + p, and deeply virtual Compton scattering (DVCS) , 

7*P ^ 7+P) have been extensively studied at HERA. These processes provide a valuable probe of the generalised (or 
00 ' skewed) gluon density at small values of the proton's momentum fraction x Possible future measurements of the 
, exclusive production of heavy vector mesons and bosons in pp collisions at the Tevatron and pp collisions at the 
' LHC should extend and complement the existing measurements from HERA in terms of both the energy span and the 
(-H , precision. In particular, by the measurement of exclusive T or J/^ production at the LHC through the detection of 
O i' two rapidity gaps with the lepton pair from the meson decay identified in the central ATLAS or CMS detectors, one 
' ' could probe the gluon density down to a; ^ 10^^ for T or x ^ 10"^ for J/ip. Even smaller values of x may be probed 
through more forward production at ALICE or LHCb. Moreover, the event rates at the LHC are expected to be 
(— I much higher than at HERA. Exclusive production at the LHC has two main merits: the process is perturbatively 
I— ' calculable with relatively small uncertainties and the experimental signature is very clean. Therefore it provides an 
interesting probe of the interplay between strong and electroweak interactions in the diffractive channel. 

The estimates of this paper are based on the equivalent-photon approximation combined with a dipole model used 
to compute the diffractive scattering of a quasireal photon on the proton. In the presence of a hard scale, such as the 
] heavy quark (or Z") mass, the impact parameter dependent dipole saturation ("b-Sat") model ^5, 's^ incorporates the 
T-H ■ impact factor for the photon to heavy meson (or transition in accordance with leading-order (LO) /cj^-factorisation, 
' together with LO DGLAP evolution of the gluon density. The input gluon density was fitted to HERA data on the 
inclusive proton structure function i^2- The model has been found to reproduce the main features of the HERA data 
for exclusive diffractive J/ip-, 4' and p production, and also for DVCS, as a function of the photon virtuality, Q^, the 
OO photon-proton centre-of-mass energy, W, and the squared momentum transfer at the proton vertex, t 0,11,11]. Thus, 
' in this paper we will use the b-Sat model to make predictions for the processes /iift.2 hi + E + h2 (E = J/ip,T, Z") 
at the Tevatron and LHC, where the reaction proceeds via photon-Pomeron fusion; see Fig.[T] The quasireal photon 
can be emitted from either of the two incoming hadrons hi. 

In Sec. [H] we recall the main formulae of the b-Sat model and in Sec. IIIII we describe the wave functions used in 
the calculation. We present the photon-proton cross sections in Sec. IIVI and the hadron-hadron cross sections in 
Sec. |Vl Finally, we compare our results with some other recent calculations in Sec. IVII and conclude in Sec. IVIII In 
the appendix we derive the light-cone wave functions for a timelike virtual photon or Z'^ boson. 

II. IMPACT PARAMETER DEPENDENT DIPOLE SATURATION MODEL 

In this section we recall the main formulae involved in the application of the b-Sat model iSj to the description 
of exclusive diffractive photoproduction. The differential cross section for the exclusive process ^p E + p is 
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FIG. 1: Exclusive photoproduction of vector mesons {J /ijj, T) or bosons in hadron-hadron interactions. The diagram where 
the photon and the Pomeron are interchanged must also be included. The variables in parentheses are the corresponding 
four-momenta. 



Here, z is the fraction of the photon's light-cone momentum carried by the quark, r — \r\ is the transverse size of 
the qq dipole, while b is the impact parameter, that is, b is the transverse distance from the centre of the proton to 
the centre of mass of the qq dipole; see Fig. [TJ The transverse momentum lost by the outgoing proton. A, is the 
Fourier conjugate variable to the impact parameter b, and t ~ — A^. We assume s-channel helicity conservation in 
the 7 ^ _E transition. We assume that the photon virtuality ^ Af|., then the contribution from longitudinally 
polarised photons may be neglected. The factor exp[i(l — z)r • A] in ([2|) originates from the nonforward wave fimctions 
0. The factor -y/Jl + 0^) 'm ^ is. d. correction to account for the imaginary part of the 5'-matrix element for the 
dipole-proton scattering, and is calculated using 

X 51n(dagg/d2b) 

The impact parameter dependent differential dipole cross section for the qq pair to scatter elastically off the proton 
is i 



d^b 



1 - exp ( -^^r'^asifi^) Rg xg{x, //^) T{b) 



(4) 



Here, the scale /i^ is related to the dipole size r hy fi^ — + /Iq. The gluon density, xg{x,iJ?), is evolved from 
a scale /Xq up to /^^ using LO DGLAP evolution without quarks. The factor Rg in @ accounts for the skewedness 
effect, that is, a; 7^ a;' in Fig. [1] and is calculated using Q 

2^^-+3 r(A. + 5/2) 9 In [xg(x,^^)] 

R^(^X,) = -^-^^—^, with A. = ^^-^^^. (5) 

The definitions of A and A^ given by ([3]) and ([5]) , respectively, are formally equivalent in the colour transparency limit 
(r — > 0), and are numerically very similar for the observables computed in this paper. ^ The initial gluon density at 
the scale /ig is taken in the form 

xg{x,^il)=AgX-^^{l~xf■''. (6) 

The values of the parameters /Iq, Ag, and Xg were determined from a fit to HERA F2 data [3]. The dipole cross 
section is evaluated at x = M'^/W'^ for E = J/t/j, T, Z'-^ The proton shape function T{b) takes a Gaussian form, that 



^ Note that the definition of A given by (|3]l diff'ers from that used in Ref . Q , where A was calculated as the logarithmic derivative of the 
7P amplitude rather than the dipole cross section. The definition of A in Jsj is more convenient for the case of photoproduction 
where there are both real and imaginary parts to the 7 — > impact factor. 
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Meson Mv/GeV /v/GeV m//GeV 


Mt 


Nl 




/y.T/GeV 


J/ip 3.097 0.274 1.4 


0.578 


0.575 


2.3 


0.307 


T(IS') 9.460 0.236 4.5 


0.469 


0.469 


0.55 


0.252 


T(1S) 9.460 0.236 4.2 


0.481 


0.480 


0.57 


0.238 



TABLE I: Parameters of the "boosted Gaussian" vector meson wave functions; see Ref. 0| for their definitions. For T the 
parameters are given for two different values of the bottom quark mass mt. 



T(6) = ^e-^, (7) 

where Bq = 4 GeV^^ is determined by the comparison to data for the t dependence of exclusive J/ if) photoproduction 
at HERA 0]. Note that although the b-Sat model incorporates saturation effects via the eikonalisation of the gluon 
density in ^ , these saturation effects are expected to be only moderate for J production and negligible for T and 
production, since the scattering amplitudes are dominated by increasingly small dipole sizes with increasing mass 
of the produced particle. 



III. PHOTON, MESON AND Z° WAVE FUNCTIONS 

The forward overlap function between the transversely polarised photon and vector meson wave functions in ([2|). 
(^I'yVE'-y)^, is given in Ref. We use the "boosted Gaussian" vector meson wave functions 0, B which were 
found to give the best description of HERA data [3] . The parameters for the J/ -0 and T wave functions are given in 
Table H 

For production, the amplitude ^ involves a sum over quark flavours f = u,d, s, c, b. The wave functions for an 
incoming Z'^ with spacelike virtuality — —Q^ < are known from the application of the colour dipole picture to 
charged-current deep- inelastic scattering [l^, [ll| . The wave functions for an outgoing Z° with timelike — > 
are derived in the appendix. The transversely polarised overlap function between the (quasireal) photon wave function 
and the wave function for quark flavour / is given by 

ino^-^yr - 4^ {[z^ + (1 - zf] mfK,imfr)ezK,iezr) + m}Ko{mfr)Ko{ezr)} . (8) 

Here, the vector couplings are g"'"^ = 1/2 — 4/3 sin^ 6w and gjj^*'^ = —1/2 + 2/3 sin^ 9w where 9w is the Weinberg 
angle, and 



Jm^r- M^z(l- z) : m'^ ~ M"^ z(l - z) > 

ez^{^ [ — — S z \ ! _ 

• > Ml z{\ ~ z) - : Mj z{l - z) - > 



The default quark masses are taken to be mu,d,s — 0.14 GeV, nic — 1.4 GeV and nib — 4.5 GeV. We use a fixed 
value of Q!cm = 1/137 for the numerical results presented in this paper, although it may be more appropriate to use a 
running acm{MD ~ 1/128, in which case the Z° cross sections would increase by 15%. 

The quasireal photon-proton cross section for production at a centre-of-mass energy W is therefore identical to 
the cross section for timelike Compton scattering, jp — s- 7*p, at a produced photon virtuality q^ — M^, apart from 
the replacement of the coupling in the amplitude: 

ee/ .^li . (10) 

Note that the timelike Compton scattering process, 7p — > 7*p, has so far only been studied at LO in the coUinear 
factorisation framework in terms of the generalised quark distribution . The timelike Compton scattering process 
at the LHC will be sensitive to the generalised gluon distribution [l3[ , and this process is calculable within the dipole 
picture using the wave functions for a timelike virtual photon given in the appendix. 

When evaluating the modified Bessel functions of an imaginary argument it is convenient to use the following 
relations of Bessel functions, valid for a real variable x > 0: 

Ko{-ix) = -^[ioW - iMx)i K,{-tx) = -|[Ji(a;) + lY.ix)]. (11) 
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FIG. 2: Choice of integration contour C to evaluate the amplitude for 7p Z^' + p when jezj > m/. 



Note that the wave function of the timelike vector boson essentially differs from the wave function of the virtual 
spacelike boson. The origin of this difference is kinematic: this point is explained in detail in the appendix. As a 
consequence, the photoproduction amplitude does not equal the electroweak DVCS amplitude at = M|. In 
particular, one sees that the overlap function, given by ([5]), picks up an imaginary part related to the contribution of 
an on-shell quark-antiquark pair at the Z'^ vertex. 

On substituting the overlap function ([5]) into the amplitude , one finds that the integrand is wildly oscillatory as 
a function of r if \ez\ ^ mf, meaning that direct numerical integration over r is difficult. This problem can be solved 
by taking the analytic continuation to complex r. By observing that the integrand is much better behaved under the 
replacement r ir one can choose the integration contour shown in Fig. [21 If there are no poles inside the contour 
C, then application of the residue theorem gives 

0=idr/(r)= / dr/(r)+ / d6l ii?e'^/(i?e'®) + / d(ir)/(ir). (12) 

Jc Jo Jo JiR 

Taking the limit R ^ oo gives the result that 

"OO />OC 

dr /(r) = / dr i/(ir). (13) 
Jo 

We use this technique to evaluate the integral over r for the case oi \ez\ > mf. (If je^j < to/, the integral over r can 
be done in the usual way. Note that there is an integrable singularity at sz — 0.) This technique requires an analytic 
form of the dipole cross section as a function of r. This can be obtained by fitting the quantity 

^Jd^te----^./irTW), (14) 

for fixed values of x and A, to a polynomial of degree 15 in log]^Q(r) for r G [10~^, 10^] GeV"'^. (Too few terms in the 
polynomial will mean that the form of the dipole cross section is not well reproduced, while too many terms will mean 
that numerical rounding errors become sizeable when r — > ir is taken due to partial cancellation between the different 
terms.) The application of equation ([13]) for a fixed z = 0.5 is illustrated in Fig. [3l It is seen that the integrand is 
much better behaved under r ^ ir such that numerical integration over r is then straightforward. 



IV. PHOTON-PROTON CROSS SECTIONS 



In Figs. [H [5] and [6| we show the jp cross sections calculated using ([TJ, integrated over \t\ < 1 GeV^, for W up to 
the maximum value of ^/s = 14 TeV theoretically accessible at the LHC. We indicate the values of the photon-proton 
centre-of-mass energy W probed by central production at the Tevatron and LHC. The predictions are compared to 
the available HERA data for exclusive J/ tp [13, dl and T [H, [13, [ll] photoproduction. For J/ tp production there is 
good agreement of the predictions with the data, as already observed in Ref . [H , but to obtain optimum agreement 
with the HERA data for the purposes of providing predictions for the Tevatron and LHC we scale the predictions for 
the "fp cross section by an overall normalisation factor of 1.08. 
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FIG. 3: Application of contour integration (I13|l to evaluate the imaginary part of the down quark contribution to the amplitude 
for 7p ^ +p for a fixed z = 0.5 with A = and x = Mzj^fs corresponding to central production at the LHC. 

For T production, however, the predictions lie roughly a factor 2-3 below the data points for our default choice 
of bottom quark mass m,\, = 4.5 GeV; see Fig. [S] The corrections for skewedness and the real part of the amplitude 
are large, at about a factor 2-3 for T, but these are already included in the "b-Sat" predictions. The T predictions 
are found to be sensitive to the assumed value of the bottom quark mass, as seen by the alternative predictions 
shown in Fig. [5] with mf, = 4.2 GeV corresponding to the running MS mass, mb{mb) |19j . Even with this low mass 
value, however, the theory curves lie significantly below the data. Such discrepancy of theory predictions and the 
HERA T data is typically found in models that use LO accuracy and an explicit meson wave function, constrained 
by the leptonic decay width; see, for example, Ref. [l^]. The successful Martin-Ryskin-Teubner (MRT) description 
of HERA T data obtained in Ref. \27\, although similar in spirit, is distinct from our approach: the && — s- T transition 
is modelled by MRT using either the meson distribution amplitude or the parton-hadron duality hypothesis; see also 
Ref. From the discrepancy between the MRT results and those from our approach, one may conclude that the 
T cross section is strongly sensitive to the choice of the T wave function. Indeed, it was demonstrated [20] that 
alternative wave functions to the "boosted Gaussian" used here give a large spread in the predictions. In addition, in 
our analysis we do not account for the higher-order QGD corrections to the impact factor, which can be quite large 
[23l[2^. All these effects, however, are expected to mainly alter the overall normalisation, but not the W dependence 
which is given by the gluon density or dipole cross section. Therefore, rather than confront these issues, we simply 
rescale the b-Sat model predictions with mi, = 4.5 GeV by a factor 2.96 to provide optimum agreement with the 
available HERA data. The rescaled prediction is also shown in Fig. [51 

As well as the b-Sat model predictions, we will also show the results of simply fitting the HERA J/tp and T data 
to the powerhke form cr(7P V + p) ex , which gives cr{jp J/iIj + p) = (3.0 nb)(W^/IVo)°'''^ and cr{jp 
T +p) = (0.12pb)(M^/M^o)^ '' with Wq = 1 GeV, shown by the dot-dashed curves in Figs. H and [5l It is clear that 
any extrapolation based on the four imprecise T data points will have a large uncertainty. In particular, note that 
the value of W probed at central rapidity at the LHC is outside the HERA kinematic range. Note also that these fits 
are to more precise HERA data than were available in Ref. [2^ . 

To investigate the numerical impact of the timelike kinematics on the cross section for Z'^ we also calculated the 
cross section with (incorrect) spacelike kinematics, that is, with — > — M| in then the amplitude is the same 
as that for DVGS at a scale = Af| apart from the different coupling. The magnitude of the imaginary part of 
the amplitude is very similar in both the timelike and spacelike cases: in the timelike case for central production it is 
0.5% smaller at the Tevatron and 2.8% larger at the LHG compared to the spacelike case. However, there is also a 
significant real part of the amplitude in the timelike case which is 24% (Tevatron) and 38% (LHC) of the imaginary 
part. The cross sections at y = are therefore enhanced by 5% at the Tevatron and 21% at the LHC in the timelike 
case compared to the spacelike case. 

For the convenience of possible future studies, we provide a simple parameterisation of the b-Sat model predictions 
for the 7p cross sections shown in Figs. HI [5] and [H] corresponding to central production at the Tevatron and LHC. 
In Table HIl we give the values of = Wq a,t y — 0, the power S ~ dlncr'^P /dlnW\w=Wo characterising the W 
dependence, and the i-slope parameter Bd obtained by fitting da'^^/dt cx exp{—B£)\t\) for \t\ < 1 GeV^. A reasonable 
approximation of the W and t dependence of the jp cross section in the vicinity of Wq may therefore be obtained 
from 



A comment on the uncertainty of this parameterisation is in order. For Jf-ip production, the accuracy of the HERA 
data and the agreement between these data and the b-Sat model suggests that the uncertainty of the theory predictions 
for the jp ^ J/4> +P cross section is 0(10%). At this level of accuracy, one should explicitly impose the rescattering 




(15) 



6 



Y p ^ J/\|/ + p 




10 10^ 10^ 10"* 

W (GeV) 

FIG. 4: W dependence of the 7p cross section for exclusive photoproduction of JM mesons. The "b-Sat" model predictions 
are rescaled by a factor 1.08 to give optimum agreement with the HERA data [13, Also shown is a direct fit to the HERA 
data of the form cr('yp) oc . The values of W corresponding to central production at the Tevatron and LHC are indicated. 
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Y p^T(1S) + p 




10 10^ 10^ lO"* 

W (GeV) 

FIG. 5: W dependence of the 7p cross section for exclusive photoproduction of T mesons. The "b-Sat" model predictions with 
nib = 4.5 GeV are rescaled by a factor 2.96 to give optimum agreement with the HERA data [l^ . [itI . [l^ . Also shown is a 
direct fit to the HERA data of the form (t(7p) oc . The values of W corresponding to central production at the Tevatron 
and LHC are indicated. 
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Y p ^ Z° + p 




10^ 10^ 10"* 

W (GeV) 

FIG. 6: W dependence of the yp cross section for exclusive photoproduction of Z° bosons in the (correct) timelike case and 
the (incorrect) spacelike case. The values of W corresponding to central production at the Tevatron and LHC are indicated. 
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Jl^, y = 


W^o (GeV) 


crTf(Wo) (nb) 


5 


Bd (GeV-2) 


Tevatron 


78 


68 


0.80 


4.66 


LHC 


208 


142 


0.71 


4.72 



T(1S), y = 


tVo (GeV) 


G^<-^{W^) (pb) 


5 


Bd (GeV-^) 


Tevatron 


136 


360 


1.39 


4.12 


LHC 


364 


1233 


1.16 


4.15 



Tevatron 
LHC 



(GeV) 
423 
1130 



°(W^o) (fb) 
4.2 
37 



3.03 
1.73 



Bu (GeV-^) 
4.25 
4.17 



TABLE II: Values of the "b-Sat" model predictions for the 7p cross section corresponding to central production at the Tevatron 
and LHC. The J j"^ and T predictions have been scaled by factors 1.08 and 2.96, respectively, in order to give the best agreement 
with the existing HERA data; see Figs. |4] and [5] A reasonable approximation of the W and t dependence in the vicinity of Wq 
may be obtained from dcr^^/dt = cr^*'(H^o) (W/Wo)''Bd exp(-Bi5ltl). 



correction of ^ 0.7-0.9 when applying the parameterisation (|15p to hadron-hadron collisions; see the later discussion 
in Sec. IVII For T production, the large experimental errors on the HERA data points and the spread between 
the various theory predictions suggests that the normalisation uncertainty factor is about 2-3; however, the energy 
dependence is expected to be accurately predicted. The estimate for production is theoretically cleanest: here we 
expect the relative uncertainty on the predictions to come mostly from higher-order QCD corrections, so to be 0(ols\ 
bearing in mind, however, that the numerical prefactor of ag in the next-to-leading-order correction may easily be 
greater than one. 



V. RAPIDITY DISTRIBUTIONS AT THE TEVATRON AND LHC 



To obtain the hadron-hadron cross sections from the photon-proton cross sections, we need to multiply by the 
photon flux dn/dfc and integrate over the photon energy k [2^: 



(j{hih2 ^ hi+E + h2) = 2 dk ^ a{-/p E + p). 

Jo dfc 



(16) 



The initial factor of 2 in accounts for the interchange of the photon emitter and the target. We neglect the 
absorptive corrections due to spectator interactions between the two hadrons and will comment on the effects of these 
in Sec. IVII The possible interference between photon-Pomeron and Pomeron-photon fusion has a large effect only 
for very small meson transverse momenta [1^ and may be safely neglected in our analysis; see also Ref. [131 • 

The four-momentum of the exchanged photon in Fig. [T] is g = (fc, qj_, fc//3i), where k and are the energy 
and transverse momentum of the quasireal photon in a given frame, where the projectile moves with velocity /3l (28| |. 
Therefore, the photon virtuality is — — ~ — fc^/(7£ (31) — q]^, where 7^ = (1 — /3|)^^/^ = ^/s/{2mp) is the Lorentz 
factor of a single beam. The photon energy spectrum is given by a modified equivalent-photon (Weizsacker- Williams) 
approximation [25l[29j: 



dn _ acm 
dfc ~ 27rfc 



2k 



In A 
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1 



A 2^2 ' 3^2") 



(17) 



where A — 1 + (0.71 GeV^)/(32^i,j and Q^jju ^ h\- This result p?]) has been obtained by integrating over the 
product of the photon propagator, l/Q^, and the squared electromagnetic form factor of the proton, F^{Q'^) 

12 //n -71 r'„^72N\-4 



+ /{0.71 GeV^))"'* [29|. We neglect the virtuality of the quasireal photon wherever possible in the calculation 
of the 7p subprocess. The square of the 7p centre-of-mass energy, 14^2 ~ 2k^/s, where ^/s is the hadron-hadron 
centre-of-mass energy. The produced state with mass Me has rapidity y ~ ln(2A:/M£;), so (fTB]) can be rewritten as 

M 



da 
dy 



{hih2 -^hi + E + h2) 



^^cr(7p 



E + p) + (y^-y), 



(18) 
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dcr/dyiy=o (nb) 


(J (nb) 


Event rate (s ^) 


Tevatron 
LHC 


3.4 
9.8 


28 
120 


0.33 
71 




T(IS') 


dcr/dyly=o (pb) 


a (pb) 


Event rate (hr ^) 


Tevatron 
LHC 


14 
72 


115 
1060 


2.0 

946 






dcr/dy|y=o (fb) 


a (fb) 


Event rate (yr ^) 


Tevatron 
LHC 


0.077 
1.4 


0.30 
13 


0.065 
135 



TABLE III: "b-Sat" model predictions for J/ip, T and photoproduction at the Tevatron Run H {y/s — 1.96 TeV) and the 
LHC (-^i = 14 TeV) using the equivalent- photon approximation. The J/^p and T predictions have been scaled by factors 
1.08 and 2.96, respectively, in order to give the best agreement with the existing HERA data; see Figs. |4]and[5l The cross 
sections must additionally be multiplied by the appropriate leptonic branching ratio for the decay E — > . These factors 
have been included when calculating the event rates, which assume a luminosity £ = 2 x 10'^'^ cm~^ at the Tevatron and 
C = 10'^'' cm~^ s~^ at the LHC. No gap survival factor has been applied to these predictions. 

where the photon energy k ~ (Mb/2) exp(y). Neglecting interference, the contribution from the interchange of the 
photon emitter and the target can be obtained by replacing y — > —y. 

In Table Hill we give the hadron-hadron cross sections at central rapidity for the Tevatron and LHC, the total cross 
sections integrated over rapidity, and the total event rates including the appropriate leptonic branching ratios and 
assuming the Tevatron and LHC design Imninosities. The cross sections for Z'^ production are comparable with the 
cross section predictions for exclusive diffractive Higgs [Mh = 120 GeV) production of 0.2 fb (Tevatron) and 3 fb 
(LHC) W{. 

In Figs. [3 [SI and [5] we show the rapidity distributions calculated using Note the sensitivity of the T rapidity 

distribution at the LHC in Fig. [5] to the W dependence of the 7p cross section: both the parameterisation and 
the rescaled b-Sat predictions describe the existing HERA data well, but give very different rapidity distributions. 
Indeed, measurements of this distribution will provide an important constraint on the generalised gluon density [sH - 

For a given rapidity y, there are contributions from proton momentum fraction x = [Me/ \/s) exp(±j/). At central 
rapidity, the x values probed at the LHC (Tevatron) are 2 x 10"^ (2 x lO"'"^) for J/i/-, 7 x 10"'' (5 x lO"-"^) for T, and 
7 X 10"'^ (5 x 10~^) for . Therefore, apart from Z° production at the Tevatron, the x values sampled for central 
rapidity are safely in the region 10"'' ^ a; < 10~^ where the dipole cross section ^ was fitted to HERA F2 data. 
However, moving further away from central rapidity, the dipole cross section must not only incorporate the correct 
dynamics at very small x, but also give sensible behaviour as a; — > 1. As noted in Ref. [25| . the simple parameterisation 
cr(7p) cx gives rise to a discontinuity at threshold; see the dot-dashed curves in Figs. [7] and [51 The b-Sat model 
is well behaved for large x, due to the damping factor (1 — a;)^'^ in the input gluon distribution Note, however, 
that this form of the input leads to a very strong x variation of the gluon density as x — > 1. Obviously, this variation 
is not related to the evolution of the gluon density in the small-x limit. Therefore, in calculating the skewing and 
real part corrections, we divide by a factor (1 — x)^'^ when calculating the value of A in ([3]) and ^ to approximately 
cancel the corresponding factor in the input gluon distribution ([5]); we set A = if the value obtained lies outside the 
interval between and 1. 



VI. DISCUSSION 

A similar approach to that used in this paper, that is, using the equivalent-photon approximation with the photon- 
proton cross sections given by the dipole model, has been applied in Ref. [s^l to calculate exclusive J/V' and T 
photoproduction at the LHC. In that case, the impact parameter dependent version of the colour glass condensate 
(CGC) model [s^ was used for the dipole cross section. This "b-CGC" model has now been updated and discussed 
further |4|]. It is less successful than the "b-Sat" model in describing observables sensitive to relatively small dipole 
sizes, such as F2 at high and J/^ photoproduction at HERA, hence we do not use it in this paper. The b-CGC 
model gives a less steep W dependence than the b-Sat model, particularly for T, and we find that we cannot reproduce 
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^ pp ^ p + J/\|/ + p at \fs = 1 .96 TeV 




pp ^ p + J/\|/ + p at \fs = 1 4 TeV 




y 



FIG. 7: Rapidity distributions for exclusive photoproduction of J/ip mesons at the Tevatron and LHC. The "b-Sat" model 
predictions are rescaled by a factor 1.08 to give optimum agreement with the HERA data [isj. Also shown is the result of 
a direct fit to the HERA data of the form (j{'yp) oc . No gap survival factor has been applied to these predictions. 
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FIG. 8: Rapidity distributions for exclusive photoproduction of T mesons at the Tevatron and LHC. The "b-Sat" model 
predictions with rUb — 4.5 GeV are rescaled by a factor 2.96 to give optimum agreement with the HERA data [l^. [itI. [l8l| . Also 
shown is the result of a direct fit to the HERA data of the form (j{'yp) oc . No gap survival factor has been applied to these 
predictions. 
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FIG. 9: Rapidity distributions for exclusive photoproduction of bosons at the Tevatron and LHC in the (correct) timelike 
case and the (incorrect) spacelike case. No gap survival factor has been applied to these predictions. 
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the results of Ref. 's?].^ Recently, an analysis of exclusive T photoproduction at hadron colliders has been performed 
within the k i -fa ctorisation framework in momentum space [2Q] . We find good agreement between our results and 
those of Ref. \M for the jp ^ T + p cross section under similar model assumptions. 

Early calculations of exclusive Z° photoproduction were made in Refs. [ssLISa]. More recently, a paper has appeared 
[sj taking a similar approach as in the present paper. The total cross sections in Ref. [sj are about a factor 5 larger 
than those presented here: such a large difference seems rather hard to explain by different model assumptions.'^ In 
that paper, another variant of the CGC model [SS^ with a <-dependent saturation scale ,39] was used. This model (and 
also the b-CGC model) are not well behaved for large x: the scattering amplitude does not vanish for a; ^ 1, and also 
the effective anomalous dimension = d\naqq/dlnr^ blows up in this limit. Hence it was necessary in Ref. (37| 
to introduce a phenomenological correction factor, (1 — x)^, where x = {Mz + /W'^. Note that Ref. [s^l did 
not include the correct timelike kinematics, but rather used the incorrect spacelike formula. However, the numerical 
impact of this mistake is only about 20% at the LHC, as shown in Fig.l^l 

Exclusive production of J /ip and T mesons via photon-Pomeron fusion was calculated within the fc_L-factorisation 
framework in Ref. [40[ as a background to odderon-Pomeron fusion. However, the main goal of Ref. [i^l was to 
provide an estimate of the oddcron contribution to exclusive diffractive vector meson hadroproduction, and not to 
provide precise predictions for the photon-mediated process. The estimates of the photon contribution were made 
there in order to monitor the impact of model assumptions in a relatively well-known situation. The "odderon to 
photon" ratio for exclusive production is expected to be strongly suppressed relative to the cases of exclusive J /'ip 
or T production due to the smaller value of the strong coupling as{M'^) compared to a5(m^) or ^^(to^), and also 
due to a stronger Sudakov suppression for glueproduced bosons because of the larger hard scale. 

It should be noted that absorptive corrections are expected to dramatically alter the ti distributions, where ti = 
{Pi — P/)^, although this effect is washed out to a large extent in the ti-integrated cross sections. Indeed, it was 
proposed in Ref. [4l| that the measurement of exclusive photon-exchange processes would provide a detailed probe 
of the rapidity gap survival probability. It was estimated in Refs. [20l . [27| that absorptive corrections reduce the ti- 
integrated rapidity distributions by a factor '-^ 0.7-0.9. Unfortunately, the inclusion of absorptive corrections requires 
knowledge of the transverse momenta of the outgoing protons, or equivalently the photon virtuality. However, this is 
already integrated over in the photon flux of (|17p . Therefore, it is not possible to include absorptive corrections, or 
to present more detailed distributions involving the momenta of the final state particles, while maintaining the simple 
equivalent-photon approximation adopted here. Instead, the rapidity distributions presented in this paper provide 
a reference point in the absence of absorptive corrections, and we refer to Refs. 20, 27, 31, 41] for a discussion of 
these effects, and for more differential distributions in the case of J/tp (and T) production. In principle, the methods 
of these papers could be applied using the b-Sat dipole model for the amplitude A{'^p — s- i? including proper 

treatment of the photon polarisation, rather than relying on an extrapolation of HERA data. 



VII. CONCLUSIONS 



We have made predictions for the rapidity distributions of exclusive photoproduced J /ip and T mesons, and Z^ 
bosons, expected at the Tevatron and LHC. We used the equivalent-photon approximation with the photon-proton 
cross sections given by the impact parameter dependent dipole saturation model [3, The normalisation of the 
J /ip and T predictions has been adjusted to give the best description of the available HERA data. For the Z^ case 
we cannot rely on existing data; however, the b-Sat model already describes well the similar process of DVCS at 
HERA [1,11]. We pointed out a crucial difference between the amplitude for DVCS and the amplitude for exclusive 
Z'^ photoproduction, and we have derived the wave functions for the case of a timelike Z'^ boson or virtual photon. 
We have given a simple parameterisation of the photon-proton cross sections for use in future studies, and we have 
discussed the uncertainties inherent in our predictions. 

Work is in progress on measurements of exclusive J /'p and T production at the Tevatron [4^ . Exclusive photopro- 
duction processes will be important at the LHC, particularly if proton taggers are installed in the 420 m region [43|. 
Exclusive production of T mesons should be measured in the early days of LHC running, extending and improving the 
existing data from HERA, and providing valuable constraints on the shape of the dipole cross section or generalised 
gluon density [3l|,[4j]. It is very unlikely that exclusive Z'^ photoproduction will be observed at the Tevatron, and it 
remains to be seen whether this process is observable at the LHC given the low event rate. 



^ The results of Ref. [s^l for T are also inconsistent with the dipole model calculations compared to HERA T data in Ref. [34| . 
^ In fact, an erroneous additional factor of 2-k in Ref. [37j has since been discovered ISSil . 
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APPENDIX: DERIVATION OF ^77* WAVE FUNCTIONS 

The dipole model approach has so far been used to describe processes with a photon or heavy boson having 
negative or vanishing virtuality, ~ —Q^ < 0. In the case of photoproduction, jp +p, or timelike Compton 

scattering, 7P ^ 7* + p, it is needed to generalise the formula to the case of = A'P > 0. 

In the dipole model the light-cone wave function of the vector boson V is defined by elementary diagrams entering 
the transition amplitude V{q)g{k) qipi)qiP2), where q, k, pi and p2 denote the particle four-momenta. The 
kinematics of this process at large energies in light-cone coordinates [for instance, p = {p^,p~,p), where =p^±p^] 
are defined as 

(7=(g+,MVg+,0), p,^(^zq+,?^^-^,p,y p,= (^(l-z)q+^^^I^,p^y (A.l) 

and k — pi + P2 — q- The three polarisation vectors of the boson are 

e± = (0,0,e±) where e±=T^(l,±i), (A.2) 

and 

where €± and ep describe the transverse and longitudinal polarisations respectively. The calculation is performed in 
the high energy limit, that is, in the limit where q^ is much larger than all other scales. 

The vector boson wave function in momentum space may be obtained from a single Feynman diagram describing 
the V{q) g{k) — > qf{pi) 9/(^2) subamplitude (see, for example, Refs. (isl. |46|) in which the gluon couples, for instance, 
to the quark line: 



[q~P2j -mj + ie 



where Pi — —P2 = P , Ca is a Dirac matrix characterising the coupling of the quark line to the incoming particle, 
and the spinors are taken in the helicity basis. The normalisation factor Af depends upon the convention: we will fix 
it later on so that the wave functions obtained match the standard dipole model expressions for the virtual vector 
boson with q^ < 0. The factor zq+ in the numerator comes from the eikonal coupling of the gluon to the quark. ^ 

Let us first analyse the case of the timelike virtual photon production with q^ — > 0. For the photon, the 
operator Oj^ takes the following form: 

61 = -iee/e^7M- (A.5) 

In the high energy limit, the matrix element (pi) O]^ f ^3 (^2) for transverse polarisations, e±, does not depend on 
the virtuality of the incoming photon, so this part of the wave function is the same for the spacelike and timelike 
virtual photons. The virtual quark propagator, however, is sensitive to the photon virtuality. In the timelike case it 
takes the form 

(g - p2)^ - m) + ie = [-p^ - m) + z{l - z)M^ + it] , (A.6) 



* In the case of the gluon coupled to the antiquark, the numerator of ||A.4| | changes into —(1 — z)q'^ , and the denominator into (g — pi) 

m) 



^ ^ ie, but the final result remains the same up to a minus sign that comes from the scattering part and not the wave function part. 
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to be compared with the standard expression for the photon with 



< 0: 



1 



{q - P2f -mj +ie = Y— -'m-} - ^{1 - z)Q^ + ie] . 



(A.7) 



In the spacehke photon case one defines a variable £^ = Q^z{l — z) + > 0. An analogous variable in the timelike 
photon case reads = rni — M^z(l — z) + ie, and its sign is not positive definite; for a large M/mf ratio will be 



negative, except in the vicinity of the z end-points. 

The wave function in coordinate space is obtained by taking the Fourier transform: 



Using the standard identities: 



d2p 



exp(ip- r)^'^i^'(p. 



(27r)2 



-ft'o(ar) for Re a >0, 



(A., 



(A.9) 



and 



d^p p fip) exp(ip • r) = -i— / d^p f{p) cxp(ip • r) 



(A.IO) 



and imposing the normalisation convention used in Refs. 0, 01 where a factor of l/(47r) appears in the integration 
measure, see we find the expressions for the wave function of the transversely polarised virtual photon: 



'i'x=±iir,z,Q) = -e/e {±ie±''l^<5Ai.±<5A..T - (1 - z)5x,.^Sx,,±]dr + mfSx,^±Sx,^±} 

for the spacelike photon, and 

^ll%{r,z,M) = -efe^/W, {±ie^'''''[zdx,^±Sx,^^ - (1 - z)dx,^^Sx,,±]dr + mfSx,^±5x,^±} 



Kojer) 
2ti 



Kojer) 
27r 



for the timelikc photon, where 



\Jmj - lVPz{l - z) : ~ APz{1- z) > 



iJAPz{l- z)-mj : M^z{l - z) - mj > 



(A.ll) 



(A.12) 



(A.13) 



and where drKn{er) = —eKi{er). Note that the signs of the various terms of (jA.lip differ from the wave functions 
given in Refs. [11,01, but the result for the overlap function summed over all helicities is the same. 

In the case of the longitudinal polarisation of the virtual photon the quark propagator is the same as above. The 
vertex part, however, is now sensitive to the incoming photon virtuality q^. One obtains 

zuA,(pi)£g7M^A2(P2) _ Vzjl^z) + mj - zjl - z)Q^ 



(q - P2f - mj +ie 



Q 



p'^ + mj + z{l - z)Q'^ 



Q 



1 - 



2z(1-z)Q2 



p2 + to2 _|_ _ 2)Q2 



(A. 14) 



for the spacelike photon, and 



z ux, jpi) s^j^ vx, iP2) _ Vzil~ z) p^ + m} + zjl - z)M^ 
[q - P2Y ~ m'j +ie 



M 



p2 + 777,2 _ — z)A'P — ie 



3Ai,-A2 



M 



2z(l - z)AP 



p'^ + mj - z{l - z)M^ ~ ie 



(A.15) 



for the timelike photon. As usual, the factors \/ z{l — z) appearing in the wave functions of the initial and final state 
are absorbed into the phase space integrations in the impact factor. The constant terms (equal to 1) in the square 
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brackets cancel in the calculation of the impact factor due to gauge invariance (the contributions of the quark and 
antiquark scattering enter the impact factor calculation with the opposite phase). It is now straightforward to obtain 
the final results for the longitudinally polarised photon: 



*aLo (r-, Q) - e/e V^<5a„-a. 2Qz[1 ~ z) (A.16) 

for the spacelike photon, and 



vf^Lo'l^. M) = -e/e ^cSm,-x, 2Mz{1 - z) (A.17) 
for the timelike photon. 

In the case of the boson, the coupling to quarks contains both vector and axial-vector parts. The amplitude for 
the Z'^ QfQf transition is described by 

r(ZO(A) ^ qfQf) = --r^e';^uf [5^7^ -5^7^75 ]«/, (A.18) 
sm Zuw 

thus 

=—^^xi9h,-9[l,l5], (A.19) 
sm Zow 

where Uf and Vf are the spinors of the quark and antiquark of the flavour /, and 6w is the Weinberg angle. The 
vector couplings are 

g-'-'t = ^-ism^ew and g,^.^." = -1 + ^ sin^ (A.20) 



while the axial-vector couplings are 



fa'''^l and 9i^^'' = -l. (A.21) 



It is natural to decompose the Z'^ wave function into distinct vector and axial-vector parts: 

vf^i^^ (r, z, M) = V^'^' (r, z, M) - A^^^' (r, z, M). (A.22) 
The analysis of the virtual photon case gives immediately the vector part of the timelike Z'^ wave function: 



sm ZUw ZTT 



sin 26 



w 



M 



2tt 



(A.24) 



The axial- vector part may be obtained in a similar way as the vector part; one needs to take into account that 
75MA = 2A?;_A (where A = ±1/2). Thus, one gets 



Ai\^^r,z,M)= - 



sin 26* 



As 



(r,z,Af) 



w 



{-ie±'^'-[z(5Ai,±<^A.,T + (1 - z)Sx,,^Sx,,±]dr ± m/(l - 2z)5x,^±Sx,^±} 



Kojir) 
2tt '■ 



(A.25) 



sin 2ew M 



{<5a„-a,2Ai [2M^z{l-z) +2m}]-iSx,,x,e-'^^''^2mfdr}^^^. (A.26) 



2tt 



The analogous formulae for the spacelike is obtained by the replacements M ^ —Q and e ^ e, then one recovers 
the formulae of Fiore and Zoller [TO, [IH up to an overall factor of 1/ y/An which in our case is included in the integration 
measure rather than the wave function; see ([2|). 

The contribution to photoproduction in the forward direction comes only from the vector part of the Z^ wave 
function. This follows from the fact that the vector and axial-vector currents have opposite parity. Therefore, the 
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overlap function for photoproduction {^*^o'^^)i^ follows from the expression for timelike Compton scattering 
i"^* '•(aP)'^i)t after the appropriate coupling adjustment, eey — s- ^.^^^ — . 
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